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Transverse energy (-Et) distributions have been measured for Au+Au collisions at ^snn = 200 
GeV by the STAR collaboration at RHIC. Et is constructed from its hadronic and electromagnetic 
components, which have been measured separately. Et production for the most central collisions is 
well described by several theoretical models whose common feature is large energy density achieved 
early in the fireball evolution. The magnitude and centrality dependence of Et per charged particle 
agrees well with measurements at lower collision energy, indicating that the growth in Et for larger 
collision energy results from the growth in particle production. The electromagnetic fraction of the 
total Et is consistent with a final state dominated by mesons and independent of centrality. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION exploration of strongly interacting matter at very high 

High energy nuclear collisions at the Relativistic Heavy 
Ion Collider (RHIC) Q have opened a new domain in the 
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energy density. High temperatures and densities may 
be generated in the most central (head-on) nuclear colli- 
sions, perhaps creating the conditions in which a phase 
of deconfmed quarks and gluons exists 0, 0] • The fireball 
produced in such collisions undergoes a complex dynam- 
ical evolution, and understanding of the conditions at 
the hot, dense early phase of the collision requires under- 
standing of the full reaction dynamics. 

Transverse energy Et is generated by the initial scat- 
tering of the partonic constituents of the incoming nuclei 
and possibly also by reinteractions among the produced 
partons and hadrons 0, @ • If the fireball of produced 
quanta breaks apart quickly without significant reinter- 
action, the observed transverse energy per unit rapidity 
cIEt I dy will be the same as that generated by the initial 
scatterings. At the other extreme, if the system rein- 
teracts strongly, achieving local equilibrium early and 
maintaining it throughout the expansion, dEx/dy will 
decrease significantly during the fireball evolution due to 
the longitudinal work performed by the hydrodynamic 
pressure 0, Q • This decrease will however be moderated 
by the buildupof transverse hydrodynamic flow, which 
increases Et |8j. Finally, gluon saturation in the wave 
function of the colliding heavy nuclei can delay the on- 
set of hydrodynamic flow, reducing the effective pressure 
and thereby also reducing the difference between initially 
generated and observed Et 0- 

Et production in nuclear collisions has been studied at 
lower at the AGS and CERN [IS [H El E H and 
at RHIC E3 • Within the framework of boost-invariant 
hydrodynamics , these measurements sug gest that en- 
ergy densities have been achieved at the SPS [Tj| that ex- 
ceed the deconfinement energy density predicted by Lat- 
tice QCD [13. However, from the foregoing discussion it 
is seen that several competing dynamical effects can con- 
tribute to the observed dEr/dy. While the measurement 
of Et alone cannot disentangle these effects, a systematic 
study of Et together with other global event properties, 
in particular charged multiplicity and mean transverse 
momentum (pt) , may impose significant constraints on 
the collision dynamics Q- 

In this paper we report the measurement of Et dis- 
tributions from Au+Au collisions at ^snn = 200 GeV 
per nucleon-nucleon pair, measured by the STAR detec- 
tor at RHIC. Et is measured using a patch of the STAR 
Electromagentic Calorimeter, with acceptance < ?/ < 1 
and A<fi = 60°, together with the STAR Time Projection 
Chamber. Et is separated into its hadronic and elec- 
tromagnetic components, with the latter dominated by 
7T° and r\ decays. The centrality dependence of Et and 
Et per charged particle are studied, and comparisons are 
made to models and to measurements at lower energy. 

A high temperature deconfmed phase could be a signif- 
icant source of low to intermediate pt photons ^3 . An 
excess of photons above those expected from hadronic 
decays has been observed at the SPS for pt > 1-5 GeV/c 
[l9j . We investigate this effect through the study of the 
electromagnetic component of Et- 



Section II describes the experimental setup used for the 
analysis. Section III presents the analysis of the hadronic 
component of the transverse energy. In section IV, the 
analysis of the electromagnetic transverse energy is pre- 
sented. In section V we discuss the scaling of Et with 
the energy of the colliding system and the number of par- 
ticipants N part and binary collisions N^n |20j, together 
with theoretical expectations for this scaling. We also 
discuss the behavior of the electromagnetic component 
of the transverse energy with the collision energy and 
centrality. Section VI is a summary and discussion of 
the main results. 



II. STAR EXPERIMENT 

This analysis is based on 150K minimum bias Au+Au 
collisions measured by the STAR detector in the 2001 
RHIC run. STAR 21] is a large acceptance, multi pur- 
pose experiment comprising several detector systems in- 
side a large solenoidal magnet. In the following, we de- 
scribe the detectors which are relevant to the present 
analysis. 

The barrel Electromagnetic Calorimeter (EMC) |22] is 
a lead-scintillator sampling electromagnetic calorimeter 
with equal volumes of lead and scintillator. It has a ra- 
dius of 2.3 m and is situated just inside the coils of the 
STAR solenoidal magnet. The electromagnetic energy 
resolution of the detector is SE/E ~ 16%/y/E(GeV). 
The results presented in this work used the first EMC 
patch installed for the 2001 RHIC run consisting of 12 
modules, ~ 10% of the full planned detector, with cov- 
erage < rj < 1 and A0 = 60°. Each EMC module 
is composed of 40 towers (20 towers in 77 by 2 towers 
in 4>) constructed to project to the center of the STAR 
detector. The transverse dimensions of a tower are ap- 
proximately 10 x 10 cm 2 , which at the radius of the front 
face of the detector correspond to a phase space interval 
of (Ary, A<^)=(0.05, 0.05). The tower depth is 21 radi- 
ation lengths (Ao), corresponding to approximately one 
hadronic interaction length. When fully installed, the 
complete barrel will consist of 120 modules with pseudo- 
rapidity coverage — 1 < rj < 1 and full azimuthal cover- 
age. 

The Time Projection Chamber (TPC) [23| has a pseu- 
dorapidity coverage of |?/| < 1.2 for collisions in the cen- 
ter of STAR, with full azimuthal coverage. In this work, 
the acceptance of the measurement was limited by the 
acceptance of the EMC. For charged tracks in the ac- 
ceptance, the TPC provides up to 45 independent spa- 
tial and specific ionization dE/dx measurements. The 
dE/dx measurement in combination with the momen- 
tum measurement determines the particle mass within 
limited kinematic regions. 

The magnetic field was 0.5 Tesla. TPC track quality 
cuts included ^-coordinate (longitudinal axis) selection of 
the collision vertex within 20 cm of the TPC center and 
a minimum TPC track space point cut of 10. Typical 
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TPC momentum resolution for the data in this work is 
characterized by 5k/ k ~ 0.0078 + 0.0098 -p T (GeV/c) || 
in which k is the track curvature, proportional to 1/pt- 
Typical resolution of dE/dx measurement is ~ 8%. Ad- 
ditional discussion of TPC analysis is given in the follow- 
ing sections and a more detailed description of the TPC 
itself can be found in Ref. |23T |. 

The event trigger consisted of the coincidence of sig- 
nals from the two Zero Degree Calorimeters (ZDC) [24j . 
located at < 2 mrad about the beam downstream of the 
first accelerator dipole magnet and sensitive to spectator 
neutrons. These calorimeters provide a minimum bias 
trigger which, after collision vertex reconstruction, corre- 
sponds to 97±3% of the geometric cross section <J^ e u t^ u - 
The events were analysed in centrality bins based on the 
charged particle multiplicity in \r]\ < 0.5. 

The procedures used in the analysis provide indepen- 
dent measurement of electromagnetic transverse energy 
and the transverse energy carried by charged hadrons. 
This latter quantity, corrected to take into account the 
contribution of the long-lived neutral hadrons, is desig- 
nated the hadronic transverse energy. The hadronic com- 
ponent of the transverse energy is obtained from mo- 
mentum analyzed tracks in the TPC while the electro- 
magnetic fraction is derived from the Electromagnetic 
Calorimeter data corrected for hadronic contamination 
using TPC tracking. In the following sections we de- 
scribe how each of these contributions were analyzed to 
obtain the total transverse energy Et measurement. 

III. HADRONIC TRANSVERSE ENERGY (E^ ad ) 
The hadronic transverse energy E T ad is defined as 



T?had 

Fjrp — 



hadrons 



E 



had 



Sint 



(i) 



where the sum runs over all hadrons produced in the 
collision, except ir° and r/. 6 is the polar angle relative to 
the beam axis and the collision vertex position. E had is 
defined for nucleons as kinetic energy, for antinucleons as 
kinetic energy plus twice the rest mass, and for all other 
particles as the total energy. E T is measured using 
charged particle tracks in the TPC via: 



E, 



had 



C 



tracks 



C 1 (ID,p)E track (ID,p) sin6 (2) 



The sum includes all tracks from the primary vertex in 
the ranges < 77 < 1 and A0 = 60°. Co is a correction 
factor defined as: 



C = 



face fpTCut fneutral 



(3) 



that includes the effective acceptance f acc — A0/27T, the 
correction fneutrai, for long lived neutral hadrons not 



measured by the TPC, and f PTC ut, for the TPC low mo- 
mentum cutoff. E tra ck{ID,p) is the energy associated 
with the particular track, either total or kinetic, as de- 
scribed above, computed from the measured momentum 
and particle identity (ID) as described below. The factor 
C\{ID,p) is defined as: 



Cx(ID,p) = f bg (p T ) 



1 



f not ID eS(p T ) 



(4) 



which includes the corrections for the uncertainty in the 
particle ID determination, f no tiD, momentum dependent 
tracking efficiency, eS(px), and momentum dependent 
backgrounds, fbgipr)- Next, we describe the corrections 
included in these two factors. 

Particle identification was carried out using the mea- 
surements of momentum and truncated mean specific 
ionization (dE/dx) in the TPC. For pt < 1 GeV/c, as- 
signment was made to the most probable particle type 
relative to the Bethe-Bloch expectation. Particles were 
assumed to be pions if (dE/dx) differed from this ex- 
pectation by more than three standard deviations, or if 
Pt > 1 GeV/c. The uncertainty in this procedure was 
gauged by calculating E T ad for px < 1 GeV/c both with 
the correct particle assignments and with all particles 
assumed to be pions. The ratio of these values for E T ad 
is applied as a correction for particles that can not be 
identified, yielding an overall correction factor to E T ad 
of fnotiD — 0.96 ± 0.02. Because this correction was 
calculated from low momentum particles it does not ac- 
count for the centrality variations in the particle ratios 
with pt > 1 GeV/c j25j. On the other hand, particles 
at pt > 1 GeV/c account for about 20% of the total 
number of particles. Taking into account the centrality- 
dependence increases in the ~p/ir and K/tt ratios at higher 
Pt generates a change in the estimated hadronic Et on 
the order of 2%, which is within the systematic error of 

fnotiD- 

Only tracks with a transverse momentum pt > 0.15 
GeV/c were accepted because the tracking efficiency 
drops rapidly below this value. GEANT [2(| detector 
simulations of HIJING [27| events demonstrate that this 
cut excludes 5% of the total E T ad . A correction f PTCU t 
for this effect is included in Co Taking all simulated 
tracks for pt < 0.15 GeV/c and calculating the energy 
assuming pions in two extreme cases, one with momen- 
tum p = and other with p — 0.15 GeV/c resulted in 
a variation of 3% in E T ad , which was assigned as the 
systematic uncertainty due to this correction. 

Since only primary charged tracks measured by the 
TPC are used in this analysis, we need to correct 
E T ad to include the contribution from long-lived neutral 
hadrons, principally n(n), K^, Kg and A(A). The cor- 
rection factor applied to the data, defined as fneutrai = 

), can be estimated u sin; 



E 



charged / /-recharged 



rpneutral 

H/rp 



measurements by STAR at 130 GeV [28J, |29j, |3fJ, |3J, 

We assume, based on HIJING simulations, that fneutrai 
does not change significantly from 130 GeV to 200 GeV. 
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We assume that the spectrum shape and yield for K\ are 
the same as for Kg. The same approximation was applied 
in the case of n(n), after subtraction of the contribution 
from A decays from the measured p yield, and the mea- 
sured STAR p/p ratio [28j. Using this procedure we ob- 
tained a value of f neutral = 0.81 ± 0.02. The uncertainty 
on this correction was estimated from the uncertainties 
in the measured STAR spectra. A cross check of these 
correction factors utilizing 200 GeV measurements j3^| 
generates variations well within the assigned systematic 
uncertainties. 

The correction fbgipr) for background, consisting of 
electrons, weak-decays and secondary tracks that are 
misidentified as primary, depends on the type of the track 
and is divided into two separate corrections. The hrst is 
for the electrons which are misidentified as hadrons. This 
correction was estimated using the shape of the electron 
spectrum obtained from HIJING and GEANT simula- 
tions and the absolute yield from STAR data in the region 
where electrons are identified with high purity using the 
TPC dE/dx measurements (essentially below 300 MeV). 
The second term is due to weak decays, which have been 
included in f neutral and therefore must be excluded from 
the primary track population to avoid double counting of 
their energy. In this case, the correction factor was calcu- 
lated by embedding simulated particles into real events. 
By comparison between the simulated particles and the 
reconstructed ones, the fraction of secondary tracks as- 
signed as primary was evaluated. A and K° were sim- 
ulated using the experimental yield and spectral shape 
measured by STAR jHHI. 

The TPC reconstruction efficiency, eff(pT)> was also 
determined by embedding simulated tracks into real 
events and comparing the simulated input and the final 
reconstructed event. In order to evaluate the effect of dif- 
ferent particle species in the reconstruction efficiency, pi- 
ons, kaons and protons were embedded in the real events. 
In this work, the charged track efficiency correction is 
the average, weighted by the relative populations of each 
of these species. The track reconstruction efficiency de- 
pends on the transverse momentum of the tracks and the 
total track density. For central events the efficiency is 
about 0.7 for tracks with px = 0.25 GeV/c and reaches a 
plateau at about 0.8 for px > 0.4 GeV/c. This efficiency 
correction includes the efficiency for track reconstruction, 
the probability for track splitting, ghost tracks and dead 
regions of the TPC. 

The resulting systematic uncertainties, taking into ac- 
count all corrections, combine in quadrature to a system- 
atic uncertainty estimate of 6.1% on E T ad . In Table[l]we 
summarize all individual corrections and the correspond- 
ing systematic uncertainties. 

HIJING and GEANT simulations of E T ad measured in 
the acceptance of this study generate event-wise fluctua- 
tions of about 10%. Simulations utilizing a substantially 
larger acceptance (0 < ?/ < 1, < <j> < 2n) generate 
event- wise fluctuations of about 4%, with this latter res- 
olution resulting mainly from tracking efficiency and neu- 



TABLE I: Corrections and systematic uncertainties for 
hadronic energy E^ ad for the 5% most central collisions. The 
quadrature sum of all the systematic uncertainties results in 
a total of 6.1%. The upper part of the table shows the global 
corrections included in Co and the bottom part shows track 
wise corrections included in Ci(ID,p). In this case, the cor- 
rection values for pr = 0.25 GeV/c and 1.0 GeV /c are shown. 



Description 


Correction 


fp T cut 


0.95 ± 0.03 


f neutral 


0.81 ± 0.02 


fnotID 


0.96 ± 0.02 


fbg(PT) 


0.84 ± 0.02 (0.25 GeV/c) 




0.94 ± 0.02 (1.0 GeV/c) 


eff(p T ) 


0.70 ± 0.04 (0.25 GeV/c) 




0.80 ± 0.04 (1.0 GeV/c) 



tral hadron corrections. 

The final E T ad distribution is corrected for vertex re- 
construction efficiency. Peripheral events have lower ver- 
tex reconstruction efficiency which suppresses the the 
transverse energy distributions with respect to more cen- 
tral events. The vertex reconstruction efficiency depends 
on the number of tracks measured in the TPC and varies 
from 70% to 97%. 



IV. ELECTROMAGNETIC TRANSVERSE 
ENERGY (E^ m ) 

The electromagnetic transverse energy E T m is the sum 
of the measured transverse energy of electrons, positrons 
and photons. The largest fraction of this energy comes 
from 7T° decays. Electrons (and positrons) are included 
because more than 90% of them are produced in the con- 
version of photons in detector materials. The energy of 
photons and electrons is fully measured by the calorime- 
ter. There is also a contribution from charged and neutral 
hadrons produced in the collision that is significant and 
must be subtracted to permit a measurement of E T m . In 
order to remove the hadronic contribution from the mea- 
surement, we studied the full spatial profiles of energy 
deposition by identified hadrons in the EMC. An exten- 
sive experimental library of hadronic shower clusters in 
the calorimeter has been obtained which, in conjunction 
with TPC tracking, allow a correction for the hadronic 
background in the calorimeter. 

Section IV. A discusses the calibration of the EMC us- 
ing minimum ionizing particles and electrons, while sec- 
tion IV. B discusses the correction for hadronic energy 
deposition in the EMC and section IV. C discusses the 
determination of E T m . 

A. Calibration of EMC 

Hadrons striking the EMC deposit a widely fluctuating 
fraction of their incident energy through hadronic show- 
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ers. In addition, ~ 30 — 40% of all high energy charged 
hadrons penetrate the entire depth of the EMC without 
hadronic interaction. If such a non-showering primary 
charged hadron has sufficient momentum, it will behave 
like a minimum ionizing particle (MIP) as it transits each 
of the scintillator layers, resulting in uniform total en- 
ergy deposition which will be nearly independent of the 
incident momentum but will vary linearly with the total 
thickness of scintillator traversed. Due to the projective 
nature of the detector the total length of the scintilla- 
tor increases with increasing r\. The MIP peak therefore 
varies from 250 MeV at small r\ to 350 MeV at large 77. 
The absolute energy of the MIP peak and its r\ depen- 
dence was determined from cosmic rays and test beam 
measurements |34j |. 

The use of MIP particles to calibrate the EMC in situ is 
convenient and provides a precision tool to track the cali- 
bration of the detector over time. In a procedure to min- 
imize systematic uncertainties in the calibration, tracks 
with p > 1.25 GeV/c in the TPC from relatively low 
multiplicity events are extrapolated to the EMC towers 
where they are required to be isolated from neighboring 
charged tracks in a 3x3 tower patch (A77 x A<f> = 0.15 
x 0.15) which has a minimum size of ~ 30 cm x 30 
cm (77 = 0). Fig. n shows a typical MIP spectrum 
measured under these conditions using minimum bias 
Au+Au events. This example shows the pseudorapidity 
interval 0.2 < 77 < 0.3. Similar spectra are observed in all 
77 bins and provide an absolute calibration in the energy 
range less than ~ 2 GeV, with an estimated systematic 
uncertainty of ~ 5% [34]]. 

An absolute calibration over a much wider energy 
range is obtained using identified electrons tracked with 
the TPC. This was done by selecting high momentum 
(1.5 < p < 5.0 GeV/c) electrons reconstructed in the 
TPC. Electron candidates are selected by dE/dx mea- 
surement in the TPC. Although the purity of the elec- 
tron candidates sample in this momentum range is poorer 
than for low momentum, the hadronic rejection factor 
obtained from the TPC dE/dx provides a clear electron 
signal in the calorimeter. Bethe-Bloch predictions for 
dE/dx of electrons and heavy particles show that the 
main background in this momentum range comes from 
deuterons and heavier particles as well as the tails of 
the distributions of protons and lower mass particles. In 
order to minimize systematic uncertainties in this proce- 
dure, only tracks having number of space points greater 
than 25 were used, as such "long tracks" exhibit better 
dE/dx resolution. It was also required that the track 
should be isolated in a 3x3 tower patch in the calorime- 
ter. 

As the final electron identifier, the energy, Etower, de- 
posited in the tower hit by the track is compared to the 
momentum, p, of the track in the range 1.5 < p < 5.0 
GeV/c. Fig. |2] shows the p/E toW er spectrum for the elec- 
tron candidates in which it is possible to see a well defined 
electron peak. The residual hadronic background in this 
figure can be evaluated by shifting the dE/dx selection 
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FIG. 1: Typical MIP spectrum. The hits correspond to iso- 
lated tracks with p > 1.25 GeV/c which project to EMC 
towers. The peak corresponds to the energy deposited by 
non-showering hadrons (MIP peak). 



window toward the pion region. The resulting estimate 
of the hadronic background is shown as a dashed line in 
the figure. After hadronic background subtraction, the 
peak position is still not centered at 1 due to the energy 
leakage to neighboring towers that is not taken into ac- 
count in this procedure. The amount of leakage depends 
on the distance to the center of the tower hit by the elec- 
tron and will shift the peak position to higher values as 
this distance increases. As shown in Fig. this effect is 
reproduced well by the full GEANT simulations of the de- 
tector response when it is hit by electrons in the momen- 
tum range used in this calibration procedure. The upper 
plot of Fig. 13 shows the position of the electron p/E tower 
peak as a function of this distance. The solid line is a 
prediction from GEANT simulations. The measurements 
are in good agreement with the simulations. Fig. flower 
plot shows the energy deposited in the calorimeter tower 
as a function of its momentum for electrons in the case 
where the distance to the center of the tower is smaller 
than 2.0 cm (A distance of 5 cm corresponds to the bor- 
der of a tower at 77 = 0. The border of the tower at 
77 = 1 is located 7.5 cm from its center). The first point 
is the electron equivalent energy of the minimum ionizing 
particles. A fit to the data using the second order poly- 
nomial of type f(x) = ao + a\x + a-ix 1 is represented by 
the solid line. The coefficients are ao = 0.01 ± 0.08 GeV, 
d = 0.98±0.11 c and a 2 = 0.01±0.03 (GeV/c 2 )' 1 . The 
values of ao and ai are consistent with zero within er- 
rors. The small magnitude of these errors indicates that 
the detector response to electrons is very linear up to 
p = 5 GeV/c. 

By combining the MIP calibration and the electron 
calibration of the EMC we obtain an overall estimated 
systematic uncertainty of less than 2% on the total en- 
ergy measured by the calorimeter. The stability of the 
detector response was evaluated by monitoring the time 
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FIG. 2: p/Etomer spectrum for electron candidates, selected 
through dE/dx from the TPC, with 1.5 < p < 5.0 GeV/c. A 
well defined electron peak is observed. The dashed line cor- 
responds to the hadronic background in the di?/da;-identificd 
electron sample. 

dependence of the shape of the raw ADC spectra for each 
tower, which is the tower response for all particles that 
reach the calorimeter. In order to have enough statistics, 
each time interval is larger than one day of data taking 
but smaller than 2 days, depending on the beam inten- 
sity during that period. The overall gain variation of the 
detector was less than 5% for the entire RHIC run. The 
results reported in this paper were obtained from 3 con- 
secutive days of data taking, in order to minimize any 
uncorrected effect due to gain variations in the detector. 

B. Energy Deposited by Hadrons in the EMC 

As discussed above, for the purposes of measuring elec- 
tromagnetic energy production it is essential to subtract 
the hadronic energy deposition in the calorimeter. For 
charged hadrons, the hit locations on the calorimeter are 
well determined and if isolated, a cluster of energy is 
readily identified. In the dense environment of Au+Au 
collisions, however, it is difficult to uniquely identify the 
energy deposition associated with a specific hadron track. 
In this limit, which is relevant for the present measure- 
ment, we subtract an average energy deposition based 
on the measured momentum of the impinging track. Be- 
cause we are interested in the cumulative distribution av- 
eraged over many events and because each event contains 
many tracks, this averaged correction results in a neg- 
ligible contribution to the uncertainty in the measured 
electromagnetic energy. 

We have studied hadron shower spatial and energy dis- 
tributions in the calorimeter both experimentally, using 
well tracked and identified hadrons in sparse events in 
STAR, and in detailed GEANT simulations. 

A library of separate profiles for pions, kaons, protons 




FIG. 3: Upper plot: points are measured p/E t0 mer electron 
peak position as a function of the distance to the center of 
the tower. The solid line is from a calculation based on a 
full GEANT simulation of the detector response to electrons. 
Lower plot: points show measured energy deposited by elec- 
trons in the tower as a function of the momentum for distances 
to the center of the tower smaller than 2.0 cm. The first point 
is the electron equivalent energy of the minimum ionizing par- 
ticles. The solid line is a second order polynomial fit of the 
data. 

and antiprotons was obtained from GEANT simulations 
of detector response in the STAR environment (GSTAR). 
The input events had a uniform momentum distribution 
in the range < p < 10 GcV/c and an emission vertex 
limited by \z verte x\ < 20 cm. The constraint on the lon- 
gitudinal coordinate of the vertex insures that the trajec- 
tory of particles will extrapolate through only one tower 
of the EMC. Because the EMC is a projective detector, 
this constraint on the extrapolated track is strongly re- 
lated to the vertex constraint. We projected the sim- 
ulated tracks on the EMC using a helix model for the 
particle trajectory in a magnetic field and obtained the 
energy distributions and the corresponding mean values 
as a function of the momentum, the pseudorapidity of the 



8 



EMC towers and the distance of the incident hit point to 
the center of the tower (d). The distributions were binned 
in intervals of A77 = 0.2. For all particles, the total mean 
deposited hadronic energy in a particular tower increases 
approximately linearly with the momentum, shows very 
little dependence on pseudorapidity and decreases with 
increasing distance from the hit point to the center of 
the tower. Experimental hadronic shower profiles were 
obtained from Au+Au minimum bias data by projecting 
tracks on the EMC, accepting only those that were iso- 
lated in a 5x5 tower patch to ensure that the energy in 
the towers was from only one particle, and calculating 
the energy distributions and mean values. Profiles for all 
particles, except electrons and positrons, for both posi- 
tive and negative tracks were recorded with good statis- 
tics up to momentum p = 2.0 GeV/c. 

In Fig. 21 we present the deposited energy for different 
particles from GEANT simulations as a function of mo- 
mentum, for a fixed pseudorapidity and distance to the 
center of the tower. An average curve, based on the rela- 
tive yield of the different particles is also presented. Small 
differences are observed for most particles, except for the 
antiproton for which the additional annihilation energy is 
apparent. The solid points are deposited energy obtained 
from experimental data for charged hadrons. The exper- 
imental profiles for charged hadrons agree quite well with 
the averaged profile. Because of the limited statistics, it 
was not possible to obtain the experimental profiles for 
identified hadrons. In Fig. [3] we present the simulated 
profiles for ir + and n~ and the experimental profiles for 
all positively and negatively charged tracks in the mo- 
mentum range 0.5 < p < 1.0 GeV/c, as a function of 
the distance to the center of a tower. The experimental 
profiles are well described by the simulation, except for a 
normalization factor on the order of 20% for < p < 0.5 
GeV/c and 5% for p > 0.5 GeV/c , as seen in Fig. 
After renormalization all experimental profiles up to mo- 
mentum p = 2.0 GeV/c are in good agreement with sim- 
ulation and we therefore use the renormalized simulated 
profiles to allow smooth interpolation in the data analy- 
sis and for extrapolation to allow corrections for higher 
momentum tracks. However, since the interval p < 2.0 
GeV/c contains 98% of all tracks, the magnitude of this 
extrapolation is small for the Et measurement. 
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FIG. 4: Mean values from GEANT simulations of the en- 
ergy deposited in the EMC by various hadronic species as a 
function of momentum. 
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FIG. 5: Spatial profiles of energy deposition in the EMC as 
a function of distance (d) from the hit point to the center of 
the tower for ir + and tv~ from simulations and for positive 
and negative hadrons from data. The arrow indicates the 
distance corresponding to the border of a tower in < rj < 0.2. 
An overall agreement between the shapes of the profiles is 
observed, with a small normalization difference (see text). 



C. E^ m measurement 



The electromagnetic transverse energy is defined as: 

EF= E EtZerM0 tower ) (5) 
towers 

where E^ er is the electromagnetic energy measured in 
an EMC tower and Qtower is the polar angle of the center 
of the tower relative to the beam axis and the collision 
vertex position. Experimentally, E^ 71 is given by: 



Ej 71 = — E {E tower — AE^ er ) Sm(8 to wer) (6) 



The sum over EMC towers corresponds to < rj < 1 
and A(f> = 60°. f acc = A0/2-7T is the correction for the 
acceptance, E toW er is the energy measured by an EMC 
tower, AE^ er is the total correction for each tower to 
exclude the contribution from hadrons . The AE!^ er 
correction is given by: 
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AE 



had 



fneutral t u C ^(Pt) 



where AE(p, 77, <i) is the energy deposited by a track pro- 
jected on an EMC tower as a function of its momentum 
p, pseudorapidity 77 and distance d to the center of the 
tower from the track hit point. feiec{pr) is a correction to 
exclude electrons that are misidentified as hadrons and, 
therefore, should not be added to AE^ er . This correc- 
tion was estimated using the same procedure described 
in the previous section to exclude real electrons from the 
E T ad measurement. eS(pr) is the track efficiency, also 
discussed previously and fneutral is the correction to ex- 
clude the long-lived neutral hadron contribution. As in 
the case for E>f d , f neutral = AE^T" l{^t^ d + 
A£ Wr a was estimated from the published STAR data 
at 130 GeV jH EE EH) In this case, AE?™%.al is 
defined as the energy deposited by all long lived neutral 
hadrons. The correction factor is fneutral = 0.86 ± 0.03. 

The systematic uncertainty due to the track effi- 
ciency correction, as previously discussed, is 4%. The 
hadronic correction for charged tracks, AE(p,r],d), is 
based primarily on measured hadronic shower profiles 
with GEANT simulations used for interpolation between 
measurements and extrapolation beyond p = 2 GeV /c. 
The systematic uncertainty for this correction to E T m 
is estimated from the observed uncertainties in the cal- 
culation of the hadronic profile at points in the shower 
library where full measurements were made. A 5% sys- 
tematic uncertainty is consistent with the comparison of 
the measured and calculated shower profiles after nor- 
malization. Different from the hadronic component of 
transverse energy, there is no correction for pt cutoff 
in the hadronic background subtraction in the electro- 
magentic energy. Such low px tracks will not reach the 
calorimeter because of the strength of the magnetic field 
and, therefore, will not deposit energy in the detector. 

As discussed earlier, the systematic uncertainty due 
to calibration of the detector is of the order of 2% and 
clearly this uncertainty contributes directly to the un- 
certainty in Elf 1 . The systematic uncertainty due to 
the electron background track correction is negligible 
(< 0.5%). 

The cumulative effect of all uncertainties discussed in 
this section, which are assumed to be uncorrelated, is 
an overall systematic uncertainty estimate for Elf of 
8.0%. All corrections and the corresponding systematic 
uncertainties are summarized in Table ITT1 

In order to evaluate the hadronic background subtrac- 
tion procedure and estimate the event-by-event resolu- 
tion of the reconstructed electromagnetic energy we have 
performed simulations in which we compare the recon- 
structed Ef 1 energy and the input from the event gen- 
erator (HIJING). Fig. (upper panel) shows the ratio, 
event by event, of the reconstructed to the input elec- 
tromagnetic transverse energy as a function of the raw 



TABLE II: Corrections and systematic uncertainties for Ef 1 
for the 5% most central collisions. The quadrature sum of all 
systematic uncertainties, including the hadronic shower pro- 
files subtraction (AE(p,rj, d)) not shown in the table, results 
in a total systematic uncertainty of 8%. The upper part of the 
table shows the global correction and the bottom part shows 
track wise corrections. In this case, the correction values for 
Pt = 0.25 Gev/c and 1.0 GeV/c are shown. 



Description 


Correction 


fneutral 


0.86 ± 0.03 


felecipT) 


0.96 ± < 0.005 (0.25 GeV/c) 




1.00 ± < 0.005 (1 GeV/c) 


eff(p T ) 


0.70 ± 0.04 (0.25 GeV/c) 




0.80 ± 0.04 (1 GeV/c) 



energy measured by the calorimeter in the same accep- 
tance used in this analysis. The smaller the raw EMC 
energy, the larger the impact parameter of the collision. 
The reconstructed energy, on average, is the same as the 
input from the event generator. Edge effects due to the 
limited acceptance of the detector were also studied and 
the effect on the reconstructed values, on average, are 
negligible. The event-by-event resolution, however, im- 
proves as the event becomes more central. Fig [B] (lower 
panel) shows the ratio distribution for the most central 
events. The solid line is a gaussian fit, from which we 
estimate that the event-by-event resolution of the recon- 
structed electromagnetic energy to be 14.5% for central 
events. The main factors that determine this resolution 
are the hadronic energy subtraction and the corrections 
for track efficiency and long-lived neutral hadrons. The 
effect on the global measurement due to the tower energy 
resolution, considering the EMC patch available, was es- 
timated to be 0.5% and that due to calibration fluctua- 
tions is 0.5%. The fluctuations due to the hadronic back- 
ground subtraction procedure alone were estimated to be 
12%, strongly dependent on the number of tracks used to 
correct the energy (for larger acceptances this resolution 
improves). The final Ef 1 distribution is also corrected 
for vertex reconstruction efficiency. 



V. TOTAL TRANSVERSE ENERGY E T 

The sum of E!f d and Ef 1 is the total transverse en- 
ergy Et of the events. In Fig. we present the Et 
distribution for minimum bias events, corrected for ver- 
tex reconstruction efficiency mainly in the low Et region. 
The scale of the upper horizontal axis corresponds to the 
Et measurement for the actual acceptance of < 77 < 1 
and A<f> — 60°. The bottom axis is scaled to correspond 
to the Et for full azimuthal coverage. 

In Fig. [7| we also present the Et distributions for dif- 
ferent centrality bins defined by the percentages of the 
total cross section, selected on charged multiplicity with 
I77I < 0.5. The centrality bin defined as 0-5% (shaded 
area in Fig. 0) corresponds to the most central collisions 
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E T m p«ch (reco) / E T . Palch (input) 



FIG. 6: Upper panel - Event-by-event ratio of the recon- 
structed electromagnetic energy and the input from the event 
generator as a function of the raw energy measured by the 
EMC. At 150 GeV, count numbers vary from 10 to 40 counts 
from the outer to the inner contour lines in steps of ~ 10 
counts. Lower panel - The same ratio distribution for the 
most central events. The solid line is a gaussian fit. 



amounting to 5% of the total cross section. The data for 
these centrality ranges are given in Table ITTT1 The cen- 
trality bins are determined by the uncorrected number 
of charged tracks with |?y| < 0.5 and number of fit points 
larger than 10. The phase space overlap between the Et 
and centrality measurements is small so that there is neg- 
ligible correlation between them beyond that due to the 
collision geometry. 

At the low energy edge, the distribution exhibits a 
peak, corresponding to the most peripheral collisions. 
For the largest values of Et the shape of the distribution 
is determined largely by statistical fluctuations and de- 
pends greatly on the experimental acceptance (35[. For 
larger acceptances, the decrease with increasing Et is 
very sharp. For this measurement, the fall off of the dis- 
tribution at large Et is strongly dominated by the lim- 
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FIG. 7: Total Transverse Energy for < 77 < 1. The minimum 
bias distribution is presented as well as the distributions for 
the different centrality bins (see table ITTTl . The shaded area 
corresponds to the 5% most central bin. The main axis scale 
corresponds to the Et measured in the detector acceptance 
and the bottom axis is corrected to represent the extrapola- 
tion to full azimuthal acceptance. 



ited acceptance which, at this point, obscures any pos- 
sible physics fluctuation. Combining the two contribu- 
tions (hadronic and electromagnetic energies) to the total 
transverse energy and properly taking into account the 
correlated uncertainties, we estimate a combined system- 
atic uncertainty in Et of 7% and an event-by-event reso- 
lution of 17%. We obtained for the 5% most central colli- 
sions (dE T /dr)\r,=o.5) = (E T ) 5 % = 621±l(stat)±43(syst) 
GeV, scaled for full azimuthal acceptance and one unit 
of pseudorapidity. 

The upper panel of Fig. [S] shows (dEx/dr]) per partic- 
ipant pair Np ar t/2 as a function of N nar t (obtained us- 
ing Monte Carlo Glauber calculations [2£j )■ Data from 
Au+Au collisions at ^/snn— 200 GeV from this analy- 
sis are shown together with similar measurements from 
Pb+Pb collisions at yi^F=17.2 GeV from WA98 [3 
and Au+Au collisions at 130 GeV from PHENIX pf . 
These comparison measurements are at r]=0, whereas the 
measurement reported here is at < r\ < 1. The grey 
bands for all three datasets show the overall systematic 
uncertainty of the data independent of N part , while the 
error bars show the quadratic sum of the statistical er- 
rors, which are typically negligible, and the systematic 
uncertainties in Et and N part |20|. with the latter dom- 
inating at low Np ar t. 

A model based on final state gluon saturation (EKRT 
7J) predicts a decrease in more central nuclear collisions 
for both the charged particle multiplicity per participant 
and Et- Hydro dynamic work during expansion may re- 
duce the observed Et relative to the initially generated 
Et, perhaps by a factor ~ 3 at RHIC energies, 7], though 
this effect will be offset somewhat by the buildup of trans- 
verse radial flow 0, @ • The dependence of observed Et 
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TABLE III: Et and E^ 1 as a function of the centrality of the collision. Global normalization uncertainties are indicated in 
the header of the table. All uncertainties are systematic. Statistical errors are negligible. 



Centrality (%) 


part 


Nbin 


Br(GeV) 
±4.3% 


£^ m (GeV) 
±4.8% 


E T /N ch (GeV) 
±5.1% 


E T /0.5N part (GeV) 
±4.3% 


rpem I jp 

Fjrp 1 11/ T 

±3.4% 


70 


-80 


14 ±4 


12 ±4 


17.1 ±0.9 


5.8 ±0.4 


0.69 ±0.07 


2.4 ±0.6 


0.342 ±0.031 


60 


- 70 


27 ± 5 


29 ± 8 




13.4 ± 0.9 


f) 75 -f- 07 


2.8 ± 0.5 


357 + 099 


50 


-60 


47 ±8 


64 ± 14 


70 ±4 


25.9 ±1.7 


0.79 ±0.06 


3.0 ±0.5 


0.369 ± 0.020 


40 


- 50 


76 ±8 


123 ± 22 


118 ±6 


43 ±3 


0.82 ±0.06 


3.1 ±0.4 


0.364 ±0.020 


30 


-40 


115 ±9 


220 ± 30 


187 ± 10 


68 ±4 


0.85 ±0.06 


3.2 ±0.3 


0.362 ±0.019 


20 


-30 


166 ±9 


368 ±41 


279 ± 15 


100 ±6 


0.86 ±0.06 


3.31 ±0.25 


0.357 ±0.019 


10 


- 20 


234 ± 8 


591 ± 52 


402 ± 21 


143 ±9 


0.86 ±0.06 


3.40 ±0.22 


0.356 ±0.019 


5- 


- 10 


299 ± 7 


828 ± 64 


515 ± 28 


181 ± 12 


0.86 ±0.06 


3.43 ±0.20 


0.351 ±0.019 





- 5 


352 ±3 


1051 ± 72 


620 ± 33 


216 ± 14 


0.86 ±0.06 


3.51 ±0.19 


0.348 ±0.019 



in y/s and system size A in the EKRT model is: 



E b T =° = 0.43A ' 92 



(^)°- 4 °(l 



0.012 \nA + 0.061 In y/s) 
(8) 

The centrality dependence can be approximated by re- 
placing A by Np ar t/2 pfl, shown by the dotted line in 
Fig El The upper panel shows comparison to measured 
cLEt I dr\ per participant pair, incorporating a Monte 
Carlo Glauber calculation for N part . The EKRT model 
is seen not to agree with the data in this panel, missing 
significantly both the centrality dependence and the nor- 
malization for central collisions. A similar comparison 
is made in the lower panel of Fig. [8] which differs from 
the upper panel only in the use of an Optical Glauber 
calculation for N part [2(]]- The centrality dependence of 
the data in this case is reproduced well by the model, 
though ~ 15% disagreement in normalization for cen- 
tral collisions remains. More precise comparison of the 
system size dependence of Erpredicted by EKRT model 
to RHIC data requires either further refinement of the 
Glauber model calculations or measurements for central 
collisions with varying mass A. 

The HIJING model predicts an increase in 
(dEx/dr]) / (0.5 Np ar t), as shown in Fig 8 upper panel. 
HIJING incorporates hard processes via the generation 
of multiple minijets together with soft production via 
string fragmentation. Effects of the nuclear geometry in 
HIJING are calculated using the Monte Carlo Glauber 
approach. Agreement of HIJING with the data is seen 
to be good in the upper panel. 

We also study a simple two-component approach where 
dEx /dij = AN part + BNt,i n . Using this model it is possi- 
ble to estimate the fraction of hard collisions in the Et 
production. In this case, a simple fit function: 



dE T /d V /(0.5N part ) = 2A[1 + (B/A){N bm /N part )} (9) 

is applied to our data at 200 GeV and the published 
PHENIX and WA98 results including points with num- 
ber of participants larger than 100. The results from the 
fits are shown in Table IIVI The simple scaling ansatz 
does a good job describing the overall shape of the N part 
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FIG. 8: (dE T /dri\r,=o.5) per N part pair vs. N part . Upper 
panel: N par t is obtained from Monte Carlo Glauber calcula- 
tions. The lines show calculations using the HIJING model 
H3 (solid) , the EKRT saturation model Q (dotted, Eq. |SJ 
and the two component fit (dashed, see text). Results from 
WA98 O and PHENIX [lj| are also shown. The grey bands 
correspond to overall systematic uncertainties, independent 
of N par t. Error bars are the quadrature sum of the errors on 
the measurements and the uncertainties on N par t calculation. 
Lower panel: the same data are shown as in the upper panel 
but using and Optical Glauber model calculation for N par t. 
The line shows the same result from EKRT model calculation. 



dependence at all energies. In this picture, the ratio 
(B/ A) (Nbin /Np a rt) estimates the fraction of the trans- 
verse energy that scales like hard processes. As seen in 
the third column in Table ITv! this ratio for the most cen- 
tral events is constant within errors despite the expecta- 
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TABLE IV: Two component model fit results of dEr/dy = 
AN pct rt ± BNbin- The uncertainties in the fit parameters in- 
clude both the data and the N par t(Nbi n ) uncertainties. 





A (Gev) 


B/A 


(B/A)(N Mn /N part ) 


STAR 


1.21 ±0.21 


0.17±0.09 


0.55 ±0.14 


PHENIX 


0.83 ±0.18 


0.27 ±0.15 


0.71 ±0.32 


WA98 


0.66 ±0.16 


0.28 ±0.11 


0.59 ±0.23 



tion that the cross section for hard processes grows by a 
large factor from 17 to 200 GeV. 

We observe an overall increase in the transverse energy 
of (24 ± 7)% at 200 GeV relative to 130 GeV. In Fig. M 
we present our result for dEr/dy per participant pair 
for central collisions, together with results from other 
experiments at various collision energies from AGS to 
RHIC 0, EE Q H3 • For the purposes of this compar- 
ison, we calculated dEr/dy from dEr/dr] for our mea- 
surements using a factor of 1.18 obtained from HIJING 
simulations to convert from r\ to y phase-space. Our re- 
sult is consistent with an overall logarithmic growth of 
dEr I ' dy I '(0.5N part ) with ^/snn- The solid line is the 
prediction using the EKRT model |7| for central Au±Au 
collisions. As one can see, EKRT model underestimates 
the final transverse energy by ~ 15%. 

We have also estimated the spatial energy density pro- 
duced in the collision using (Er)$% reported above, con- 
verted from pseudorapidity to rapidity density using the 
factor of 1.18 discussed above. Based on a scaling solu- 
tion to the relativistic hydrodynamic equations, Bjorken 
[lil estimated the spatial energy density of the system in 
terms of the primordial transverse energy rapidity density 
dEr/dy, the transverse system size, R, and a formation 
time To. 



dE 7 



' B i dy tqitR 2 



(10) 



We assumed tq = 1 fm/c which is the usual value 
taken in many analysis at SPS energies. For Au±Au 
200 GeV we obtained 



at y/sNM 



-Bj 



4.9 ± 0.3 



GeV/fm 3 . The uncertainty includes only the uncertainty 
on (dEr/dr}). This energy density is significantly in ex- 
cess of the energy density ~ 1 GeV/fm 3 predicted by 
Lattice QCD for the transition to a deconfmed Quark 
Gluon Plasma 01- The estimate is based however upon 
the assumption that local equilibrium has been achieved 
at r ~ 1 fm/c and that the system then expands hydro- 
dynamically. Comparison of other RHIC data, in partic- 
ular elliptic flow, to hydrodynamic calculations [42l |43| 
indicates that this picture may indeed be valid. 

In order to understand the systematic growth in 
transverse energy with collision energy shown in 
Fig. we investigate the centrality dependence of 
(dEr/drj)/(dN c h/drj), the scaling of transverse energy 
relative to the number of charged particles produced in 
the collision. The centrality dependence of this ratio may 
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FIG. 9: dEr/dy (see text for details) per N var t pair vs ^/snn 
for central events. In this figure dEr /dy / (0.5 N par t) is seen to 
grow logarithmicaly with ^/sjvjv- The error bar in the STAR 
point represents the total systematic uncertainty. The solid 
line is a EKRT model prediction corrected for dy/dy, for 
central Au±Au collisions. 



indicate effects of hydrodynamic flow y| : if the expansion 
is isentropic then dN c h/dr) will remain constant, whereas 
dEr /df] will decrease due to the performance of longitu- 
dinal work. 

Fig. upper panel, shows the centrality depen- 

dence of (dEr/dr]) / (dN c h/drj) from STAR measurements 
at ^snn = 200 GeV, compared to similar measurements 
at 17 and 130 GeV. Data at all energies fall on a common 
curve within uncertainties, with modest increase from 
the most peripheral collisions to -/V part =100, reaching a 
roughly constant value of (dEr/dr/)/ {dNch/drj). Fig. 1101 
lower panel, shows the (pr) for 200 GeV Au±Au col- 
lisions measured by STAR |37| . showing a dependence 
on centrality similar to that of the transverse energy per 
charged particle: modest increase with N part for N part < 
100, with constant value for more central collisions. The 
systematic behavior of Er, multiplicity and (pr) are sim- 
ilar, indicating that the growth of Er is due to increased 
particle production. Quantitative comparison of theo- 
retical models of particle production with the measured 
centrality dependences of (dEr / drf) / (dN c h/ drf) and (pr) 
of charged particles will constrain the profile of initial 
energy deposition and the role of hydrodynamic work 
during the expansion. 

In Fig. El we show, for central collisions, that this con- 
stant transverse energy production per charged particle 
is observed down to and including AGS measurements at 
^snn = 5 GeV. A single value of ~ 800 MeV per charged 
particle or at most a slow logarithmic increase amounting 
to < 10% characterizes all measurements within errors 
over a range in which the Er per participant grows by 
a factor of four. HIJING predicts that Er per charged 
particle should increase from SPS to RHIC energies due 
to the enhancement of minijet production at RHIC. How- 
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FIG. 10: Upper panel - (dE T /drf) / (dN ch /drj) vs N part . Pre- 
dictions from HIJING simulations for Au+Au at 200 GeV are 
presented. Results from WA98 [3 and PHENIX O are also 
shown. The grey band corresponds to an overall normaliza- 
tion uncertainty for the STAR measurement. Bottom panel 
- Charged hadrons mean transverse momentum as a function 
of N part Ht 



ever, the predicted increase is rather small and the sys- 
tematic uncertainties on the measurement do not provide 
enough precision to significantly test this assumption. 

The procedures adopted in this analysis permit an 
independent measurement of the electromagnetic and 
hadronic transverse energy. This allows additional explo- 
ration of the collision dynamics and particle production. 
In Fig. ^]we show the ratio of the electromagentic to 
the total energy for the most central events as a function 
of the energy from lower SPS energies |llll38l | to our re- 
sults at full RHIC energy. The observed electromagnetic 
fraction of the total transverse energy will be strongly 
influenced by the baryon to meson ratio. At very high 
energy it is expected that virtually all the Et will be car- 
ried by mesons and the fraction should approximate 1/3, 
whereas at low energy, baryon dominance of the trans- 
verse energy will result in a much smaller electromagnetic 
fraction. 

While the energy dependence seen in Fig. ^] is pre- 
sumably dominated by the total meson content of the 
final state, the centrality dependence may provide addi- 
tional detail about the reaction mechanisms. The cen- 
trality dependence of the electromagnetic fraction of our 
total measured energy is shown in Fig. 1131 An excess 
photon yield may result from the formation of a long- 
lived deconfined phase, as suggested in Ref. |39j. The 
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atic uncertaintiy. A constant value of ~ 800 MeV per charged 
particle, within errors, characterizes transverse energy pro- 
duction over this full energy range. 



predictions from HIJING simulations are also presented. 
We observe no significant dependence of the electromag- 
netic fraction with the collision centrality. 



VI. SUMMARY 

We have reported the measurement of transverse en- 
ergy Et within < r\ < 1, for centrality-selected Au+Au 
collisions at ^/snn = 200 GeV. For the 5% most central 
events we measured (Et)$% = 621 ± l(stat) ± 43(syst) 
GeV, corresponding to an increase of (24±7)% with re- 
spect to measurements at 130 GeV Au+Au collisions at 
RHIC HI. 

We investigated the energy scaling with the number of 
participant nucleons and with the number of charged par- 
ticles produced in the collision. We obtained, for the 5% 
most central events, <1Et / 'drj / '(0.5N part ) = 3.51 ± 0.24 
GeV and (dE T /dr))/{dN ch /dri) = 860 + 70 MeV, re- 
spectively. We also compared the results of this work 
with measurements from AGS and SPS energies. It was 
found that the increase in the Et production from AGS 
up to RHIC energies comes mostly from the increase in 
the particle production. A final state gluon saturation 
model (EKRT), HIJING and a simple two component 
(hard/soft) model were compared to the data. Although 
the EKRT model predicts a different centrality behavior 
of energy production the uncertainties in the N part de- 
termination does not allow us to discard this model. The 
simple two component ansatz suggests that, despite the 
large uncertainties, the fraction of energy arising from 
hard processes which is still visible in the final state does 
not increase significantly from SPS to RHIC energies. 

Other measurements at RHIC and comparison to the- 
oretical calculations suggest that a dense, equilibrated 
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FIG. 12: Energy dependence of the electromagnetic fraction 
of the transverse energy for a number of systems spanning 
SPS to RHIC energy for central events. 
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FIG. 13: Participant number dependence of the electromag- 
netic fraction of the total transverse energy. The results are 
consistent with HIJING within errors over the full centrality 
range. 



system has been generated in the collision and that it 
expands as an ideal hydrodynamic fluid. The good agree- 
ment between hydrodynamic calculations and measure- 
ments of particle-identified inclusive spectra and elliptic 



flow [42j are consistent with the onset of hydrodynamic 
evolution at a time To < 1 fm/c after the collision [43| . 
The strong suppression phenomena observed for high 
hadrons |40l 144 . |45| suggest that the system early in its 
evolution is extremely dense. Estimates based on these 
measurements yield an initial energy density in the vicin- 
ity of 50-100 times cold nuclear matter density. Within 
the framework of boost-invariant scaling hydrodynamics 
|l6j . from the Et measurement presented here we esti- 
mate an initial energy density of about 5 GeV/fm 3 . This 
should be understood as a lower bound , due to the 
strong reduction in the observed relative to the initially 
produced Ex from longitudinal hydrodynamic work dur- 
ing the expansion. These three quite different approaches 
produce rough agreement for the estimated initial energy 
density, with a value well in excess of that predicted by 
Latice QCD for the deconfinement phase transition [17| . 

The method used in this analysis permitted an in- 
dependent measurement of the electromagnetic and 
hadronic components of the total energy. The elec- 
tromagnetic fraction of the transverse energy for the 
5% most central events obtained in this work is 
(dE^/dr]) I (dE T /dr/) = 0.35 ± 0.02, consistent with a 
final state dominated by mesons. Some models [39l | ex- 
pect that the formation of a long lived deconfined phase 
in central events may increase the yield of direct pho- 
ton production and, therefore, an increase in the electro- 
magnetic fraction of the transverse energy. We, however, 
observe that the electromagnetic fraction of the trans- 
verse energy is constant, within errors, as a function of 
centrality. Measurements with larger acceptances would 
have systematic uncertainties significantly reduced and 
therefore would be able to show smaller effects that can 
not be observed with the precision of the present mea- 
surement. 
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